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Abstract. We present VLT eclipse photometry for the giant planet CoRoT-lb. We observed a transit in the 
R-band filter and an occultation in a narrow filter centered on 2.09 /im. Our analysis of this new photometry 
^ . and published radial velocities, in combination with stellar-evolutionary modeling, leads to a planetary mass and 

H ■ radius of 1.07to]_lMjup and 1.45^^;^ Rjup, confirming the very low density previously deduced from CoRoT 

C/3 , photometry. The large occultation depth that we measure at 2.09 fj.ia (0.278lg Qgg%) is consistent with thermal 



0.066 ' 

emission and is better reproduced by an atmospheric model with no redistribution of the absorbed stellar flux to 
the night side of the planet. 

Key words, binaries: eclipsing - planetary systems - stars: individual: CoRoT-1 - infrared: stars - techniques: 
photometric - techniques: radial velocities - techniques: spectroscopic 



1. Introduction Spitzer spectral window were performed (Richardson et 

^ al. 2003a,b; Snellen 2005; Deming et al. 2007; Knutson 

O Transiting planets play an important role in the study of ^1. 2007, Snellen & Covino 2007, Winn et al. 2008). 

g planetary objects outside our solar system. Not only can yery recently, two of them were successful: Sing & Lopez- 

V: ; we infer their density and use it to constraint their com- Morales (2009) obtained a - 4 a detection of the occulta- 

_> . position, but several other interesting measurements are ^-^^ OGLE-TR-56b in the z-band (0.9 , while De 

possible for these objects (see e.g. review by Charbonneau ^oiij & Snellen (2009) detected at - 6 a the thermal emis- 

. et al. 2007). In particular, their thermal emission can be ^^^^ TrES-3b in the K-band (2.2 It is important 

Ci ; measured during their occultation, allowing the study of ^^^^-^ ^^^^ gj^^jl^^, measurements to improve our un- 

their atmosphere without spatially resolving their light derstanding of the atmospheric properties of short-period 

from that of the host star. The Spitzer Space Telescope extrasolar planets 
(Werner et al. 2004) has produced a flurry of such plan- 

etary emission measurements, all at wavelengths longer CoRoT-lb (Barge et al. 2008, hereafter B08) was the 

than 3.5 ^m. From the ground, several attempts to obtain f^'^s* planet detected by the CoRoT space transit sur- 

occultation measurements at shorter wavelengths than ^^y (Baghn et al. 2006). With an orbital period of 1.5 

days, this Jupiter-mass planet orbits at only ^ 5 stel- 
lar radii from its GOV host star. Due to this proximity. 
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* T, 1 1 . „ . , . , .1 -<rrm ^r^T,r,r, i its stcUar irradiatiou IS clcarly largc cuough (~ 3.9 X IQ-^ 

* Based on data collected with the VLT/FORS2 and , ., . . rr. r^r. o, 

VLT/HAWK-I instruments at ESO Paranal Observatory, Chile ^ ) "^^'^^ JO^'^ OGLE-TR-56b, TrES-3b 

(programs 080.C-0661(B) and 382.C-0642(A)). and a few other planets withm the pM theoretical class 

** The photometric time-series used in this work are proposed by Fortney et al. (2008). Under this theory, pM 

only available in electronic form at the CDS via anony- planets receive a stellar flux large enough to have high- 

mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via opacity compounds like TiO and VO present in their 

|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-|-A/ gaseous form in the day-side atmosphere. These com- 
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pounds should be responsible for a stratospheric thermal 
inversion, with re-emission on a very short time-scale of 
a large fraction of the incoming stellar flux, resulting in 
a poor efficiency of the heat distribution from the day- 
side to the night-side and to enhanced infrared plane- 
tary fluxes at orbital phases close to the occultation. Like 
the other pM planets, CoRoT-lb is thus a good target 
for near-infrared occultation measurements. Furthermore, 
CoRoT-lb belongs to the subgroup of the planets with 
a radius larger than predicted by basic models of irra- 
diated planets (e.g. Burrows et al. 2007, Fortney et al. 
2007). Tidal heating has been proposed by several au- 
thors (e.g. Bodenheimer et al. 2001, Jackson et al. 2008b) 
as a possible extra source of energy able to explain the 
radius anomaly shown by these hyper-bloated planets. As 
shown by Jackson et al. (2008b) and Ibgui & Burrows 
(2009), even a tiny orbital eccentricity is able to pro- 
duce an intense tidal heating for very short period plan- 
ets. Occultation photometry does not only allow to mea- 
sure the planetary thermal emission, but also constrains 
strongly the orbital eccentricity (see e.g. Charbonneau et 
al. 2005). Such an occultation measurement for CoRoT-lb 
could thus help understanding its low density. 

These reasons motivated us to measure an occultation 
of CoRoT-lb with the Very Large Telescope (VLT). We 
also decided to obtain a precise VLT transit light curve 
for this planet to better constrain its orbital elements. 
Furthermore, CoRoT transit photometry presented in B08 
is exquisite, but it is important to obtain an indepen- 
dent measurement of similar quality to check its reliability 
and to assess the presence of any systematic effect in the 
CoRoT photometry. 

We present in Section 2 our new VLT data and their 
reduction. Section 3 presents our analysis of the resulting 
photometry and our determination of the system param- 
eters. Our results are discussed in Section 4, before giving 
our conclusion in Section 5. 



2. Observations 

2.1. VLT/F0RS2 transit photometry 

A transit of CoRoT-lb was observed on 2008 February 28 
with the F0RS2 camera (Appenzeller et al. 1998) installed 
at the VLT/UTl (Antu). F0RS2 camera has a mosaic of 
two 2k X 4k MIT CCDs and is optimized for observations 
in the red with a very low level of fringes. It was used 
several times in the past to obtain high precision transit 
photometry (e.g. Gillon et al. 2007a, 2008). The high res- 
olution mode was used to optimize the spatial sampling, 
resulting in a 4.6' x 4.6' field of view with a pixel scale of 
0.063"/pixel. Airmass increased from 1.08 to 1.77 during 
the run which lasted from lhl6 to 4h30 UT. The quality of 
the night was photometric. Due to scheduling constraints, 
only a small amount of observations were performed before 
and after the transit, and the total out-of-transit (OOT) 
part of the run is only ~ 50 minutes. 
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Fig. 1. Top: VLT/F0RS2 R-band transit light curve 
with the best-fitting transit + trend model superimposed. 
Bottom: residuals of the fit. 



114 images were acquired in the R_SPECIAL filter 
(Ae// — 655 nm, FWHM = 165 nm) with an exposure 
time of 15 s. After a standard pre-reduction, the stellar 
fluxes were extracted for all the images with the IRAF0 
DAOPHOT aperture photometry software (Stetson, 1987). 
We noticed that CoRoT-1 was saturated in 11 images 
because of seeing and transparency variations, and we 
rejected these images from our analysis. Several sets of 
reduction parameters were tested, and we kept the one 
giving the most precise photometry for the stars of simi- 
lar brightness than CoRoT-1. After a careful selection of 
reference stars, differential photometry was obtained. A 
linear fit for magnitude vs airmass was performed to cor- 
rect the photometry for differential reddening using the 
OOT data. The corresponding fluxes were then normal- 
ized using the OOT part of the photometry. The resulting 
transit light curve is shown in Fig. 1. After subtraction of 
the best- fit model (see next section), the obtained resid- 
uals show a rms of '~ 520 ppm, very close to the photon 
noise limit (~ 450 ppm). 



2.2. VLT/HAWK-I occultation photometry 

We observed an occultation of CoRoT-lb with HAWK- 
I {High Acuity Wide-field K-band Imager, Pirard et al. 
2004, Casali et al. 2006), a cryogenic near-IR imager re- 
cently installed at the VLT/UT4 (Yepun). HAWK-I pro- 
vides a relatively large field of view of 7.5' x 7.5'. The 
detector is kept at 75 K and is composed of a mosaic of 
four Hawaii-2RG 2048x2048 pixels chips. The pixel scale 
is 0.106" /pixel, providing a good spatial sampling even for 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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the excellent seeing conditions of Paranal (seeing down to 
0.3 arcsec measured in K-band). 

Instead of using a broad band K or filter, we choose 
to observe with the narrow band filter NB2090 (central 
wavelength = 2.095 /im, width = 0.020 /zm). This filter 
avoids absorption bands at the edge of the K-band, its 
small width minimizes the effect of differential extinction 
and furthermore its bandpass shows a much smaller sky 
emission than the one of the nearby Br7 bandpass (central 
wavelength = 2.165 /Ltm, width = 0.030 /itm), leading to a 
flux ratio background/star more than twice better than in 
Br7 or K-band filters. Because of the large aperture of the 
VLT and the relative brightness of CoRoT-1, the expected 
stellar count in this narrow filter is still good enough to 
allow theoretical noise of less than 0.15% for a 1 minute 
integration. 

Observations took place on 2009 January 06 from 
lh54 to 7h56 UT. Atmospheric conditions were very good, 
while the mean seeing measured on the images was 0.47" . 
Airmass decreased from 1.36 to 1.08 then raised to 1.65. 
Each exposure was composed of 4 integrations of 11 s each. 
A random jitter pattern within a square 45" -sized box was 
applied to the telescope. This strategy aimed to obtain for 
each image an accurate sky map from the neighbor images. 
Indeed, the near-IR background shows a large spatial vari- 
ability at different scales and an accurate subtraction of 
this complex background is crucial, except when this back- 
ground has an amplitude negligible when compared to the 
stellar count (see e.g. Alonso et al. 2008). In total, 318 im- 
ages were obtained during the run. 

After a standard pre-reduction (dark subtraction, flat- 
field division), a sky map was constructed and removed 
for each image using a median-filtered set of the ten ad- 
jacent images. The resulting sky-subtracted images were 
aligned and then compared on a per-pixel basis to the 
median of the 10 adjacent images in order to detect any 
spurious values due, e.g., to a cosmic hit or a pixel dam- 
age. The concerned pixels had their value replaced by the 
one obtained by linear interpolation using the 10 adjacent 
images. 

Two different methods were tested to extract the stel- 
lar fluxes. Aperture photometry was obtained using the 
IRAF DADPHDT software and compared to deconvolution 
photometry obtained with the algorithm DECPHOT (Gillon 
et al. 2006 and 2007b, Magain et al. 2007). Wc obtained 
a significantly (-- 25%) better result with DECPHOT. We 
attribute this improvement to the fact that DECPHOT op- 
timizes the separation of the stellar flux from the back- 
ground contribution, while aperture photometry simply 
sums the counts within an aperture. 

In order to avoid any systematic noise due to the dif- 
ferent characteristics of the HAWK-I chips, we choose to 
use only reference stars located in the same chip than our 
target to obtain the differential photometry. As CoRoT-1 
lies in a dense field of the Galactic plane, we have in any 
case enough reference fiux in one single chip to reach the 
desired photometric precision. After a careful selection of 
the reference stars, the obtained differential curve shows 
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Fig. 2. Top: VLT/HAWK-I 2.09 ^m occultation hght 
curve binned per 10 minutes, with the best-fitting occul- 
tation + trend model superimposed. Bottom: residuals of 
the fit. 



clearly an eclipse with the expected duration and timing 
(Fig. 2). We could not find any firm correlation of the 
OCT photometric values with the airmass or time, so we 
simply normalized the fluxes using the OOT part without 
any further correction. The OOT rms is 0.32 %, much 
larger than the mean theoretical error: 0.13 %. This dif- 
ference implies the existence of an extra source of noise 
of ~ 0.3 %. We attribute this noise to the sensitivity and 
cosmetic inhomogeneity of the detector combined with our 
jitter strategy. In the optical, one can avoid this noise by 
staring at the same exact position during the whole run, 
i.e. by keeping the stars on the same pixels. In the near-IR, 
dithering is needed to remove properly the large, complex 
and variable background. This background varies in time 
at frequencies similar to the one of the transit, so any poor 
background removal is able to bring correlated noise in the 
resulting photometry. It is thus much preferable to opti- 
mize the background subtraction by using a fast random 
jitter pattern even if this brings an extra noise, because 
this latter is dominated by frequencies much larger than 
the one of the searched signal and is thus unable to pro- 
duce a fake detection or modify the eclipse shape. 

3. Analysis 

3.1. Data and model 

To obtain an independent determination of the system pa- 
rameters, we decided to use only our VLT R-band transit 
and 2.09 fim occultation photometry in addition to the 
SOPHIE (Bouchy et al. 2006) radial velocities (RV) pre- 
sented in B08 as data for our analysis. 

These data were used as input into a Markov Chain 
Monte Carlo (MCMC; see e.g. Tegmark 2004, Gregory 
2005, Ford 2006) code. MCMC is a Bayesian inference 
method based on stochastic simulations and provides the 
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a posteriori probability distribution of adjusted param- 
eters for a given model. Here the model is based on a 
star and a transiting planet on a Kcplcrian orbit about 
their center of mass. More specifically, we used a classi- 
cal Keplerian model for the RV variations and we fitted 
independent offsets for the two epochs of the SOPHIE 
observations to account for the drift between them men- 
tioned in BOS. To fit the VLT photometry, we used the 
photometric eclipse model of Mandcl & Agol (2002) mul- 
tiplied by a trend model. In order to obtain reliable er- 
ror bars for our fitted parameters, it is indeed prefer- 
able to consider the possible presence of a low-amplitude 
time-dependant systematic in our photometry due, e.g. to 
an imperfect differential extinction correction or a low- 
amplitude low-frequency stellar variability. We choose to 
model this trend as a second-order time polynomial func- 
tion for both F0RS2 and HAWK-I photometry. 

3.2. Limb-darkening 

For the transit, a quadratic limb darkening law was as- 
sumed, with initial coefficients ui and U2 interpolated from 
Claret's tables (2000; 2004) for the R-band photometric 
filter and for T^fj = 5950 ± 150 K, log g = 4.25 ± 0.30 
and [Fe/H] = -0.30 ± 0.25 (BOS). We used the partial 
derivatives of ui and U2 as a function of the spectroscopic 
parameters in Claret's tables to obtain their errors (T„j 
and (T,., via: 



\ 



1=1 



(1) 



where a; is 1 or 2, while Si and as- are the i*'' (i = 1,3) 
spectroscopic parameter and its error from BOS. We ob- 
tained ui = 0.279 ± 0.033 and U2 = 0.351 ± 0.016 as ini- 
tial values. We allowed ui and U2 to float in our MCMC 
analysis, using as jump parameters not these coefficients 
themselves but the combinations ci = 2 x ui + U2 and 
C2 = Ml — 2 X U2 to minimize the correlation of the ob- 
tained uncertainties (Holman et al. 2006). The following 
Bayesian prior on ci and C2 was added to our merit func- 
tion: 



limb— darkening 



E 



(2) 



where c- is the initial value deduced for the coefficient c, 
and ac'. is its error computed from cr^j and 0-^2 • We let ci 
and C2 be free parameters imdcr the control of a Bayesian 
prior to propagate the uncertainty on the limb-darkening 
to the deduced transit parameters. 

3.3. Jump parameters 

The other jump parameters in our MCMC simulation 
were: the transit timing (time of minimum light) Tq, 
the planet/star area ratio (Rp/Rs)"^, the transit width 



(from first to last contact) W, the impact parameter 
b' = ocosi/i?*, three coefficients per photometric time- 
scries for the low-frequency systematic, one systemic RV 
for each of the two SOPHIE epochs, and the two param- 
eters ecosw and esinw where e is the orbital eccentricity 
and uj is the argument of pcriastron. The RV orbital semi- 
amplitude K was not used as jump parameter, but instead 
we used the following parameter: 



K2 = KVl-e^ = (27rG)i 



/3. 



siu i, 



(3) 



to minimize the correlation with the other jump parame- 
ters. 

We notice that our used jump parameter b' is equal to 

the actual transit impact parameter b only for a circular 
orbit. For a non-zero eccentricity, it is related to the actual 
impact parameter b via: 



6 = 6' 



1 



1 -|- esinw 



(4) 



Here too, the goal of using b' instead of b is to minimize 
the correlation between the jump parameters. 

The orbital period P was let free in our analysis, 
constrained not only with the data presented above but 
also with the timings determined independently by Bean 
(2009) for each of the 35 CoRoT transits. Practically, 
we added the following bayesian penalty BPtimings to our 
merit function: 



BPt 



timings 



E 

i=l,35 



To + N,xP-T, 



(5) 



where Ti is the transit timing determined by Bean (2009) 
for the CoRoT transit, ctt; is its error and Ni is its 
differential epoch compared to our VLT transit. This pro- 
cedure relies on the reasonable assumption that the tim- 
ings determined by Bean (2009) are uncorrelated with the 
other transit parameters. 



3.4. Photometric correlated noise and RV jitter noise 

Our analysis was done in 4 steps. First, a single MCMC 
chain was performed. This chain was composed of 10® 
steps, the first 20% of each chain being considered as its 
burn-in phase and discarded. The best-fitting model foimd 
in the first chain was used to estimate the level of corre- 
lated noise in each photometric time-series and a jitter 
noise in the RV time series. For both photometric time- 
series, the red noise was estimated as described in Gillon 
et al. (2006), by comparing the rms of the unbinned and 
binned residuals. We used a bin size corresponding to a 
duration of 20 minutes, similar to the timescale of the 
ingress/egress of the transit. The obtained results were 
compatible with a purely Gaussian noise for both time- 
series. Still, it is possible that a low-amplitude correlated 
noise damaging only the eclipse part had been 'swallowed' 
by our best-fitting model, so we preferred to be conser- 
vative and to add quadratically a red noise of 100 ppm 
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to the theoretical uncertainties of each photometric time- 
series. The deduced RV jitter noise was high: 23 m.s~^. 
Nevertheless, we noticed that it goes down to zero if we 
discard the second RV measurement of the first SOPHIE 
epoch. Furthermore, this measurement has a significantly 
larger error bar than the others, we thus decided to con- 
sider it as doubtful and to do not use it in our analysis. 
A theoretical jitter noise of 3.5 m.s~^ was then added 
quadratically to the error bars of the other SOPHIE mea- 
surements, a typical value for a quiet solar- type star like 
CoRoT-1 (Wright 2005). 

3.5. Determination of the stellar density 

Then, 10 new MCMC chains were performed using the up- 
dated measurement error bars. These 10 chains were then 
combined, using the Gelman and Rubin statistic (Gelman 
& Rubin 1992) to verify that they were well converged 
and mixed enough, then the best-fitting values and error 
bars for each parameter were obtained from their distri- 
bution. The goal of this MCMC run was to provide us 
with an improved estimation of the stellar density p* (see 
e.g. Torres 2007). The stellar density that we obtained was 
= 0.841°: pQ. 

3.6. Stellar-evolutionary modeling 

The deduced stellar density and the spectroscopic param- 
eters were then used to better constrain the stellar mass 
and age via a comparison with theoretical stellar evolution 
models. Two independent stellar analysis were performed 
in order to assess the impact of the stellar evolution mod- 
els used on the final system parameters: 

— Our first analysis was based on Girardi's evolution 
models (Girardi et al. 2000). We first perform a linear 
interpolation between the solar (Z=0.019) and subso- 
lar (Z=0.008) metallicity theoretical models to derive 
a set of mass tracks at the metallicity of the host star 
([M/H]=-0.3). We then compare the effective temper- 
ature and the inverse cube root of the stellar density 
to the same values in the host star metallicity mod- 
els. We interpolate linearly along the mass tracks to 
generate an equal number of age points between the 
zero age main sequence and the point corresponding 
to core hydrogen exhaustion. We then interpolate be- 
tween the tracks along equivalent evolutionary points 
to find the mass, M = 0.94 M©, and age, r = 7.1 
Gyr, of the host star that best match the measured 
temperature and stellar density. We repeat the above 
prescription using the extreme values of the observed 
parameters to determine the uncertainties on the de- 
rived mass and age. The large errors on the spectro- 
scopic parameters, particularly the ±0.25 dex uncer- 
tainty on the metallicity, lead to a 15-20 % error on 
the stellar mass (M = O.94lo;JgM0) and an age for 
the system no more precise than older than 0.5 Gyr. 
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Fig. 3 presents the deduced position of CoRoT-1 in a 
i?/MV3-Te// diagram. 
— In the second analysis we apply the Levenberg- 
Marquard miniminization algorithm to derive the fun- 
damental parameters of the host star. The merit func- 
tion is defined by: 

3 (r>obs /^theo\2 
y2 = ^^i - L (6) 

i=l ^ i ' 

The observables (Of^) we take into consideration are 
effective temperature, surface metallicity and mean 
density. The corresponding observational errors are 
^obs -pjjg theoretical values {Of^^°) are obtained from 
stellar evolution models computed with the code CLES 
(Code Liegois d'Evolution Stellaire, Scuflaire et al. 
2008). Several fittings have been performed, in all 
of them we use the mixing-length theory (MLT) of 
convection (Bohm-Vitense, 1958) and the most re- 
cent equation of state from OPAL (OPAL05, Rogers 
& Nayfonov, 2002). Opacity tables are those from 
OPAL (Iglesias & Rogers, 1996) for two different so- 
lar mixtures, the standard one from Grevesse & Noels 
(1993, GN93) and the recently revised solar mixture 
from Asplund, Grevesse & Sauval (2005, AGS05). 
In the first case {Z/X)q=Q.Q2A^, in the second one 
{Z / X)q=O.Q\&J . These tables are extended at low 
temperatures with Ferguson et al. (2005) opacity val- 
ues for the corresponding metal mixtures. The surface 
boundary conditions are given by grey atmospheres 
with an Eddington law. Microscopic diffusion (Thoul 
et al. 1994) is included in stellar model computation. 
The parameters of the stellar model are mass, ini- 
tial hydrogen (Xi) and metal {Zi) mass fractions, age, 
and the parameters of convection (aMLT and the over- 
shooting parameter). Since we have only three obser- 
vational constraints, we decide to fix the aMLT and Xi 
values to those derived from the solar calibration for 
the same input physics. Furthermore, given the low 
mass we expect for the host star, all the models are 
computed without overshooting. The values of stel- 
lar mass and age obtained for the two different solar 
mixtures are: M = 0.90 ± 0.21 Mq with GN93 and 
M = 0.92 ± 0.18Mq with AGS05, and respectively 
T = 7.5 ± 6.0 Gyr and r = 6.9 ± 5.4 Gyr. 

The result of our two independent stellar analysis are 
thus fully compatible, and the uncertainty due to the large 
errors on the spectroscopic parameters dominates the one 
coming from our imperfect knowledge of stellar physics. 
The large uncertainties affecting the stellar mass and age 
are mainly due to the lack of accuracy in metallicity de- 
termination. We estimate from several tests that an im- 
provement in the atmospheric parameters determination 
leading to an error in metallicity of 0.05 dex would trans- 
late in a reduction in uncertainty by a factor three for the 
stellar mass, and a factor two for the stellar age. Moreover, 
decreasing the effective temperature error to 75 K would 
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Fig. 3. R/M^/^ (in solar units) versus effective tempera- 
ture for CoRoT-1 compared to tlie theoretical stellar stel- 
lar evolutionary models of Girardi et al. (2000) interpo- 
lated at -0.3 metallicity. The labeled mass tracks are for 
0.8, 0.9 and 1.0 Mq and the isochrones are 100 Myr (sohd), 
5 Gyr (dotted), 10 Gyr (dashed), 16 Gyr (dot-dashed). We 
have interpolated the tracks at -0.2 metallicity and have 
included the uncertainty on the metallicity (±0.25) in the 
overall uncertainties on the mass and the age. 

imply a subsequent reduction of stellar parameter errors 
by an additional factor two. Getting more high-SNR high- 
resolution spectroscopic data for the host star is thus very 
desirable. 

3.7. Determination of tlie system parameters 

For the last part of our analysis, we decided to use 0.93 ± 
0.18 Mq, i.e. the average of the values obtained with the 
two different evolution models, as our starting value for 
the stellar mass. A new MCMC run was then performed. 
This run was identical to the first one, with the exception 
that M^, was also a jump parameter under the control of 
a Bayesian penalty based on Af* = 0.93 ± 0.18 Mq. At 
each step of the chains, the physical parameters Mp, Rp 
and i?* were computed from the relevant jump parame- 
ters including the stellar mass. Table 1 shows the deduced 
values for the jump + physical parameters and compares 
them to the values presented in BOS. It also shows the 
Bayesian penalties used in this second MCMC run. 

4. Discussion 

4.1. The density and eccentricity of CoRoT-lb 

As can be seen in Table 1, the transit parameters that we 
obtain from our VLT/FORS-2 R-band photometry agree 
well with the ones presented in B08 and based on CoRoT 
photometry. Our value for the transit impact parameter is 
in good agreement with the one obtained by BOS, and has 
a similar uncertainty. The planet/star area ratio that we 
deduce is within the error bar of the values obtained by 



BOS, while our error bar is smaller. Our deduced physical 
parameters agree also very well with the ones presented 
in BOS. Our analysis confirms thus the very low-density 
of the planet (see Fig. 4) and its membership to the sub- 
group of short period planets too large for current models 
of irradiated planets (Burrows et al. 2007; Fortney et al. 
2008). 

In this context, it is worth noticing the marginal non- 
zero eccentricity that we deduce from our combined analy- 
sis: e = 0.071^0 028- -^^ outlined by recent works (Jackson 
et al. 2008b, Ibgui & Burrows 2009), tidal heating could 
play a major role in the energy budget of very short pe- 
riod planets and help explaining the very low density of 
some of them. Better constraining the orbital eccentric- 
ity of CoRoT-lb by obtaining more radial velocity mea- 
surements and occultation photometry is thus desirable. 
To test the amplitude of the constraint brought by the 
occultation on the orbital eccentricity, we made an anal- 
ysis similar to the one presented in Sec. 3 but discard- 
ing the HAWK-I photometry. We obtained similar re- 
sults for the transit parameters, but the eccentricity was 
poorly constrained; we obtained much less precise val- 
ues for ecosoj and esincj, respectively 0.0201q q29 ^^.d 
-0.170lo:o78. HAWK-I occultation brings thus a strong 
constraint on these parameters, especially on ecosw. 

Table 1 shows that our analysis does not agree with 
BOS for one important parameter: the stellar density. 
Indeed, the value presented in BOS is significantly smaller 
and more precise than ours. Still, BOS assumed a zero ec- 
centricity in their analysis, while the stellar mean density 
deduced from transit observablcs depend on e and lu (see 
e.g. Winn 2009). To test the influence of the zero eccen- 
tricity assumption on the deduced stellar density, we made 
a new MCMC analysis assuming e = 0. We obtained this 
time ~ 0.695^0 030 PQ^ excellent agreement with the 
value = 0.698 ±0.033 p© presented by BOS. This shows 
nicely that not only are VLT and CoRoT data fully com- 
patible, but also that assuming a zero eccentricity can lead 
to an unreliable stellar density value and uncertainty. In 
our case, this has no significant impact on the deduced 
physical parameters because the large errors that we have 
on the stellar effective temperatm-c and metallicity dom- 
inate totally the result of the stellar-evolutionary model- 
ing (see Sec. 3.6). Still, the point is important. As shown 
by Jackson et al. (2008a), most published estimates of 
planetary tidal circularization timescales used inappropri- 
ate assumptions that led to unreliable values, and most 
close-in planets could probably keep a tiny but non-zero 
eccentricity during a major part of their lifetime. In this 
context, very precise transit photometry like the CoRoT 
one is not enough to reach the highest accurary on the 
physical parameters of the system, a precise determina- 
tion of e and co is also needed. This strengthens the inter- 
est of getting complementary occultation photometry in 
addition to high-precision radial velocities to improve the 
characterization of transiting planets. 
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Parameter 


Value 


Bayesian penalty 


BOS 


Unit 


Jump parameters 



Transit epoch To 

Planet /star area ratio {Rp/RsY 

Transit width W 

2.09 ^m occultation depth 

6' = acosi/Rt 

RV K2 

RV71 

RV 72 

ecosw 

esinw 

-'^transit 
-^transit 

a ransit 

-^occultation 
-^occultation 
^'occultation 



2454524.623241°;°°°?^ 
01<30fi+° °°°^° 



FO. 00043 
^-0.00066 
j+0.032 

-0.043 

-16 
-.+0.020 
-0.017 



0.10439 ± 0.00094 
0.002781 
0.3981 
215:1 
23.3661 

23.3501°:°!? 
0.0083l°;°°g 
-0.070l°g;°f2 

U.333UO_g ggggg 

017+0 003 

"•'-'J^' -0.018 
'-'•^'J-0.02 



1.00041 



-0.00096 
-0.00052 



-0.008l°g;°°^ 
0.029l°g:°S 



2454159.4532 ± 0.0001 
0.01927 ± 0.00058 



0.420 ± 0.043 
216 ± 13 



BJD 

day 

R* 

km s~^ 
km s~^ 



day" 
day 

day 
day 



Orbital period P 
Stellar mass M* 
R-filter ci 
R-filter C2 



1.50896861 



1 01+013 

J-'Ul.g 22 

704+0 047 

-0 444+0 054 



from timings in Bean (2009) 
0.93 ± 0.18 
0.909 ± 0.067 
-0.423 ± 0.046 



1.5089557 ± 0.0000064 
0.95 ± 0.15 



day 
M© 



Deduced parameters 



RV K 


188 ± 14 


188 ± 11 


ms ^ 


^transit 


4 26+0 035 

U.'±ZO_Q g42 


0.420 ± 0.043 




^occultation 


Q7n+0 037 

U.O(U_Q g4g 


0.420 ± 0.043 




Orbital semi-major axis a 


0259+0-00" 


0.0254 ± 0.0014 


AU 


Orbital inclination i 


85.66lg:f3 


85.1 ± 0.5 


degree 


Orbital eccentricity e 


o.o7il°:°g 


(fixed) 




Argument of periastron uj 


276.71^:;^ 




degree 


Stellar radius R* 


1 r.t;7+0.055 

l.UO(_g g94 


1.11 ± 0.05 


Rq 


Stellar density p* 


0.86t°:J^8 


0.698 ± 0.033 


Po 


R-filtcr Ml 


229+0 025 
"•^^■'-0.022 






R-filter U2 


S36+0 012 

U.OOU_Q g2g 






Planet radius Rp 


1 45+0.07 


1.49 ± 0.08 


Rj 


Planet mass Mp 


1 07+013 

■■■•"-"-0.18 


1.03 ± 0.12 


Mj 


Planet density pp 


350+0 077 

U.OOU_g g42 


0.31 ± 0.06 


PJ 



Table 1. CoRoT-1 system parameters and 1-a error limits derived from the MCMC analysis. The parameters A, B 
and C are the zero-, first- and second-order coefiicients of the polynomial used to model the photometric trend. The 
vahics and error bars used in the Bayesian penalties are shown in the third column. Fourth column shows the values 
presented in BOS. 



4.2. The atmospheric properties of CoRoT-lb 

The flux at 2.09 /im of this planet is slightly larger than 
the (zero-albedo) equilibrium temperature, ~ 2660 K, ob- 
tained if the star's effective temperature is allowed to be 



as high as 6100K (maximum within the l-cr error-bars 
from BOS). An irradiated planet atmosphere model (fol- 
lowing Barman et al. 2005) for CoRoT-lb, was computed 
adopting the maximum observational allowed stellar ef- 
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12 3 4 



M (Mj) 

Fig. 4. Position of CoRoT-lb (in red) among the 
other transiting planets (black circles, values from 
|http://exoplanet.eu) in a mass-radius diagram. The error 
bars are shown only for CoRoT-lb for the sake of clarity. 




2 3 4 5 

Wavelength {jiim) 



Fig. 5. Comparison of our 2.09 /im occultation depth mea- 
sured for CoRoT-1 with models of planet-star flux den- 
sity ratios assuming that the absorbed stellar flux is re- 
distributed across the dayside only (top curve) and uni- 
formly redistributed across the entire planetary atmo- 
sphere (lower curve). A black body model is also shown 
(dotted) for T = 2365 K. 



fective temperature and radius and assuming zero energy 
is transported to the night side. Solar metallicity was as- 
sumed and all other parameters were taken from Table 1. 
This model (Fig. 5) falls short of matching the observa- 
tions within 1-cr, while a black body with the same equi- 
librium temperature as the irradiated planet model is in 
better agreement. The atmosphere model is hot enough 
for a significant temperature inversion to form for P < 
0.1 bar and is nearly isothermal from 0.1 down to ~ 100 
bar. A model which uniformly redistributes the absorbed 
stellar flux across the entire planet surface (lower curve 
in Fig. 5) is far too cool to match the observations and is 
excluded at 3 cr. The flux at 2.09 /im alone is suggestive 
that very little energy is redistributed to the night side; 
however additional observations at shorter and/or longer 
wavelengths are needed to better estimate the bolometric 
flux of the planet's day side. Occultation measurements 
in other bands will help provide limits on the day side 
bolometric flux and determine the depth of any possible 
temperature inversion and the extent of the isothermal 
zone. 

Recently Snellen et al. (2009) measured the dayside 
planet-star flux ratio of CoRoT-1 in the optical (~ 0.7 
/im) to be 1.26±0.33 x lO"''. The hot, dayside only, model 
shown in Fig. 5 predicts a value of 1.29 ± 0.33 x 10"**, 
which is fully consistent with the optical measurement. 
Consequently, it appears as though very little energy is 
being carried over to the night side of this planet. 



4.3. Assessing the presence of another body in the 
system 

As shown in Table 1, our deduced systemic RV for each 
SOPHIE epoch agree with well each other: we do not con- 
firm the RV drift mentioned in BOS. Our combined analy- 
sis presented in Sec. 3 leads to a very precise determination 
of the orbital period: 1.5089686t°:o°o°o°g days, thanks to a 
lever arm of nearly one year between CoRoT transits and 
the VLT one. A simple linear fit for timing vs epoch based 
on the CoRoT and VLT transits lead to a similar level of 
precision, giving P = 1.5089686l|^;;^^|^J5^^^ days. This fit 
has a reduced of 1.28 and the rms of its residuals (see 
Fig. 6) is 36 s. These values are fully consistent with those 
reported by Bean (2009) for CoRoT data alone. We also 
notice the same 3-cr discrepancy with transit #23 which, 
once removed, results in a reduced of 1.00, hereby con- 
firming the remarkable periodicity of the transit signal. 

While limits on additional planetary companions in 
CoRoT-1 system have been extensively discussed for tran- 
sit timing variations (TTVs) by Bean (2009), it is inter- 
esting to take limits from RVs into account. Figure 7 il- 
lustrates the domains where additional planets could be 
found through TTVs (white) and through RV measure- 
ments (above coloured curves). We focused on short pe- 
riods objects since TTVs are more sensitive to nearby 
perturbators as compared to the known transiting planet. 
We assumed an eccentricity of 0.05 for a putative copla- 
nar planet and used the Mercury package described in 
Chambers (1999) to estimate by numerical integration the 
maximum TTV signal expected for CoRoT-lb. White is 
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Fig. 6. Top: Residuals of the hnear fit timing vs epoch 
for CoRoT-lb (see text for details). The rightmost point 
is our VLT/F0RS2 timing. Botttom: zoom on the CoRoT 
residuals. 

the domain with a > 5-cr detection through TTVs ac- 
cording to CoRoT data rms, while black area is below 
the l-fj detection threshold. Although approximative, this 
shows that for a typical 3 ms^^ accuracy on radial ve- 
locities (dashed curve on Fig. 7) routinely obtained with 
HARPS spectrograph (Mayor et al. 2003), no room re- 
mains for planetary companion detection through TTVs 
alone. Thus, TTV search method may opportunely be 
applied on active and/or stars for which RV measure- 
ments accuracy is limited, increasing a detectability area 
for which RVs arc not or far less sensitive to. 

Each transit timing may be compared to a single RV 
measurement. The increased free parameters in a TTV or- 
bital solution raise degeneracies that could not be lifted by 
considering the same number of datapoints that would al- 
low an orbital solution recovery with RVs. Determination 
of a large number of consecutive transits added to occul- 
tation timings allow to access to the uniqueness of the 
solution as well as lowering constraints on timings accu- 
racy (Nesvorny & MorbidcUi 2008). This is thus a high 
cost approach that is the most potentially rewarding on 
carefully determined targets stars. 

5. Conclusion 

We have obtained new high-precision transit photometry 
for the planet CoRoT-lb. Our deduced system parame- 
ters are in very good agreement with the ones presented 
in BOS, providing thus an independent verification of the 
validity of the CoRoT photometry. Due to the precision 
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t; 
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to 




10.0 



Period [d] 

Fig. 7. Detectivity domain for a putative CoRoT-lc 
planet, assuming Cc — 0.05. In white, the period-mass 
region where planets yield maximum TTV on CoRoT- 
lb above 100 s (Scr detection based on CoRoT data). 
Companions in the black area yield maximum TTV below 
the la threshold. Solid, dashed and dotted curves shows 
RV detection limits for 1, 3 and 10 m/s RMS. 

of the CoRoT and VLT transit photometry and the long 
baseline between them, the orbital period is now known to 
a precision better than 1/lOth of a second. The precision 
on the planetary mass and radius is limited by the large 
errors on the stellar spectroscopic parameters, and a sig- 
nificant precision improvement should be made possible 
by getting new high-quality spectra of CoRoT-1. 

We have also measured successfully the occultation of 
the planet with HAWK-I, a new wide-field near-infrared 
imager mounted recently on the VLT. The large occulta- 
tion depth that we measure is better reproduced by an 
atmospheric model with no redistribution of the absorbed 
stellar flux to the night side of the planet. This measure- 
ment firmly establishes the potential of the HAWK-I in- 
strument for the study of exoplanetary atmospheres. At 
the time of writing, Spitzer cryogen is nearly depleted 
and soon only its 3.6 /im and 4.5 fiui will remain available 
for occultation measurements, while the eagerly awaited 
JWST (Gardner et al. 2006) is not scheduled for launch 
before 2013. It is thus reassuring to note that ground- 
based near-infrared photometry is now able to perform 
precise planetary occultation measurements, bringing new 
independent constraints on the orbital eccentricity and on 
the atmospheric physics and composition of highly irradi- 
ated extrasolar planets. 
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